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Abstract 

Over the summer of 2007, we obtained 1191 hours of 'drift-scan' pulsar search obser- 
vations with the Green Bank Telescope at a radio frequency of 350 MHz. Here we describe 
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the survey setup, search procedure, and the discovery and follow-up timing of thirteen pul- 
sars. Among the new discoveries, one (PSR J1623— 0841) was discovered only through its 
single pulses, two (PSRs J1327— 0755 and J1737— 0814) are millisecond pulsars, and an- 
other (PSR J2222-0137) is a mildly recycled pulsar. PSR J1327-0755 is a 2.7 ms pulsar 
at a DM of 27.9 pc cm~^ in a 8.7 day orbit with a minimum companion mass of 0.22 M0. 
PSR J1737— 0814 is a 4.2 ms pulsar at a DM of 55.3 pc cm~^ in a 79.3 day orbit with a 
minimum companion mass of 0.06 Mq. PSR J2222— 0137 is a 32.8 ms pulsar at a very low 
DM of 3.27 pc cm~^ in a 2.4 day orbit with a minimum companion mass of 1.11 Mq. It is 
most likely a white dwarf-neutron star system or an unusual low-eccentricity double neutron 
star system. Ten other pulsars discovered in this survey are reported in the companion paper 
Lynch et al. 2012. 

Key Words: starsrneutron - pulsarsrgeneral - pulsars: surveys - pulsars:individual 
(PSR J1327-0755, PSR J1623-0841, PSR J1737-0814, PSR J1941-0121, PSR J2222-0137) 



1. Introduction 

Since the discovery of the first pulsar in 1967 (Hewish et al. 1968), many fascinating 
sources and experiments have been associated with pulsars. These include but are not limited 
to tests of general relativity (e.g. Lyne et al. 2004, Kramer et al. 2006), studies of binary 
evolution (e.g. Bhattacharya & van den Heuvel 1991), probes of electron density models of 
the interstellar medium (e.g. Cordes & Lazio 2002), constraints on the neutron star equation 
of state (e.g. Lattimer & Prakash 2001, Lattimcr & Prakash 2007), and use in a pulsar timing 
array for gravitational wave detection (e.g. Hobbs et al. 2010, Demorest et al. 2012). 

Searches for pulsars are very demanding, requiring both significant amounts of compu- 
tational and human resources. The high time and frequency resolution needed for pulsar 
searches result in data rates greater than for most astronomical projects (typically 100 GB/hr 
or more). Large computing clusters are needed to thoroughly search for both periodic and 
transient dispersed signals. After all the computer searches are completed, candidates are 
typically visually inspected to determine likely pulsar candidates, though automated algo- 
rithms are becoming more widely used for sifting through large numbers of candidate signals 
(e.g., Eatough et al. 2010). 

Pulsars are brighter at low frequencies (i.e. a few hundred MHz) but pulsar survey 
sensitivity is reduced by dispersion smearing, higher sky temperature, and scattering at these 
frequencies. Therefore most recent large-scale pulsar surveys like the Parkes Multibeam 
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Pulsar Survey (PMPS; Manchester et al. 2001), High Time Resolution Universe Pulsar 
Survey (HTRU; Keith et al. 2010), and the Pulsar Arecibo L-band Feed Array Survey 
(PALFA; Cordes et al. 2006) are conducted around 1.4 GHz. Lower frequency observations 
(< 500 MHz) are better used to locate weak sources away from the Galactic plane where the 
effects of scattering, higher sky temperature, and dispersion broadening are minimized. 

The National Radio Astronomy Observatory (NRAO) in Green Bank, WV has a long 
history of successful searches for radio pulsars at low frequencies. Damashek et al. (1978) 
discovered 17 new pulsars with the 92-meter (300') transit telescope at a frequency of 400 
MHz. Later, using the same telescope at a frequency of 390 MHz, Dewey et al. (1985) 
found 34 new long-period pulsars. However, they were not sensitive to newly discovered 
millisecond pulsars (MSPs; P < 20 ms; Backer et al. 1982). The first low-frequency survey 
to be sensitive to MSPs was conducted by Stokes et al. (1985) at 390 MHz. While this 
survey did not discover any MSPs, it did find 20 new long-period pulsars. Sayer et al. (1997) 
used the 43-meter (140') telescope at Green Bank at a frequency of 370 MHz and discovered 
eight new pulsars including one mildly relativistic binary pulsar. New backend hardware 
and higher data rates now allow frequency resolution that was previously unattainable by 
earlier observing systems. Most recently, a North Galactic plane survey was conducted with 
the Robert C. Byrd Green Bank TelescopeEl (GBT) by Hessels et al. (2008) at 350 MHz and 
discovered 33 new pulsars. None of these were MSPs despite the survey's sensitivity to these 
system^. 

Our survey took advantage of a unique opportunity when the GBT was immobilized for 
track refurbishing during the northern summer of 2007. During this time, we recorded data 
at a radio frequency of 350 MHz as the sky passed through the beam of the telescope. We 
will discuss these drift-scan observations in Section 2. In Section 3 we describe the follow-up 
observations and timing analysis of the new pulsars. In Section 4 we present the properties 
of the thirteen pulsars and in Section 5 we discuss interesting individual objects. Finally in 
Section 6 we present the conclusions and the future survey planned in light of the results 
obtained in this survey. The survey sensitivity, pipeline, and the remaining new pulsars are 
discussed in Lynch et al. (2012; Paper II). 



^The Robert C. Byrd Green Bank Telescope (GBT) is operated by the National Radio Astronomy Ob- 
servatory. 

^The North Galactic plane survey has only been partially processed to look for MSPs. 
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2. Survey Observations 

The drift-scan observations occurred from May through August in 2007. The telescope 
was parked at a set azimuth during the track refurbishing and data were recorded while the 
sky drifted through the telescope's beam. The elevation of the telescope was adjusted and 
had a limit of > 60° during the weekdays and > 30° during the weekends. The azimuth 
was set at ~229° for the first half of the observations and ~192° for the second half of 
the observations. These azimuth restrictions allowed us to access declination (6) ranges of 
-7.72° <6 < 38.43° and -20.66° <6 < 38.43° for azimuths of 229° and 192°, respectively 
All right ascensions were available for us to observe. 

We chose a frequency of 350 MHz for this survey for many reasons. Firstly, at 350 
MHz, the radio frequency interference (RFI) environment at the GET is well known and 
remarkably manageable. Secondly, pulsars generally have very steep spectral indices (e. g. 
Sieber 1973) and our survey covers Galactic latitudes away from the Galactic plane, where 
propagation effects in the interstellar medium are generally less severe. Lastly, the GBT's 
beam size at 350 MHz is ~36', allowing us to cover a greater amount of the sky than would 
be possible with higher frequency receivers. The sky coverage rate (210 deg^ per day) is 
faster than that achieved by the PMPS at 1.4 GHz (26 deg^ per day; Manchester et al. 
2001). 

The now-retired Spigot autocorrelation spectrometer (Kaplan et al. 2005) was used for 
the GET drift-scan observations. We used 50 MHz of bandwidth with 81.92-/xs sampling 
time. For the first week of observations (MJDs 54222 - 54230) data were recorded with 
1024 three-level auto-correlation functions which were accumulated into 16-bit registers for 
analysis. After that time, and until the conclusion of the observations, a new 2048 three-level 
auto-correlation mode was used which accumulated data into 8-bit registers. 

Data were recorded every night on each weekday and for 24 hours a day during the 
weekends. In total, 1491 hours of data were taken, amounting to 134 TE (see Supplemental 
Online Material for the observation list and Figured] for sky coverage area). This amounts 
to a sky coverage of 10347 square degrees, of which ~2800 went to the Pulsar Search Col- 
laboratory (PSC; Rosen et al. 2010; Rosen et al. 2012)§ A source passing through the 
center of the beam is visible for roughly 140 s. Assuming a detection significance of 8a, a 
digitization correction factor of 1.16 to account for 3-level quantization of the signal, and a 
pulse duty cycle of 0.04 results in a survey flux density limit 



■^The PSC is a collaboration between NRAO and WVU that involves high-school students in searching 
for pulsars, http: / / www.pulsarsearchcollaboratory.com/ 
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S^in = 0.012T,y,— ^, (1) 
is. 

where Tgys is the total system temperature. The system temperature has two major com- 
ponents: the receiver temperature, which is ~22 K, and the sky temperature. The sky 
temperature can vary greatly in the survey region. Far off the plane, the sky temperature is 
~30 K at 350 MHz but near the plane the sky temperature can be 200—400 K in our survey 
regions. For a total system temperature of 52 K, a survey sensitivity limit of 0.6 mJy is ob- 
tained. Scaling our flux density limit with an average pulsar spectral index of —1.8 (Maron 
et al. 2000) gives us a limit of ~0.05 mJy at 1352 MHz, more sensitive than the mid- and 
high-latitude HTRU survey in the same regions, but our observations were hampered more 
significantly by dispersion smearing and scattering. A more detailed discussion of the survey 
sensitivity is presented in Paper H. Data processing was done with the PRESTC0 software 
package (Ransom 2001) and the data pipeline will be presented in Paper H. 



3. Timing Observations and Analysis 

Timing observations began in June of 2008. At that stage, nine new pulsars that had 
been discovered were observed. The first session was conducted at a central frequency of 
350 MHz with the GET using the Spigot backend with 50 MHz of bandwidth and 81.92-/is 
sampling time. With observations spread across six hours, this session was used to establish 
a very accurate period and for gridding of pulsars to find a more accurate position (Morris 
et al. 2002). 

Most of the rest of the follow-up timing observations in the first year (July 2008 — May 
2009) were conducted using the Spigot backend at a central frequency of 820 MHz with 1024 
channels, 50 MHz of bandwidth and 81.92 /is sampling time. All further observations (June 
2009 onwards) were conducted using the Green Bank Ultimate Pulsar Processing Instrument 
(GUPPI) backend (Ransom et al. 2009). Most of the observations used a center frequency 
of 820 MHz with 200 MHz of bandwidth, 2048 channels, and 40.96 /is sampling time. For 
pulsars that had satisfactory timing solutions, the data were taken in on-line folding mode 
with ten-second integrations, while the rest were done in search mode, with the data recorded 
in PSRFITS format (Hotan et al. 2004). 

Each search-mode data set was first subbanded, or de-dispersed at the pulsar's dispersion 
measure (DM) into broader frequency channels, using PRESTO. We used 16 subbands for long- 



'*http://www.cv.nrao.edu/^sransoni/presto/ 
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period pulsars and 32 subbands for MSPs. The DM value used for subbanding is either the 
discovery value or a timing- derived value, if multi-frequency times of arrival (TOAs) had 
allowed fitting the DM. The on-line folded GUPPI data require no subbanding. 

Folded pulsar profiles were created from the subbanded data using PRESTO; an example 
output plot is shown in Figure |2l These folds were done using the best available ephemeris 
with 128 bins for normal period pulsars or 64 bins for MSPs. On-line folded GUPPI data 
were co-added in time and frequency and flux calibrated (described later) to form a profile 
using the PSRCHIVE0 package with 256 pulse phase bins. 

TOAs were measured from the folded profiles using the frequency-domain algorithm in 
PRESTO for all search-mode data (Taylor 1992). We measure one TOA per observation for 
each isolated pulsar and three per observation for each binary pulsar for better determination 
of binary parameters. TOAs for the on-line folded data were created with PSRCHIVE and also 
used the Taylor (1992) method with the same number of TOAs as the search-mode data. 
Templates were created by fitting multiple Gaussian functions to the summed pulse profile 
for each frequency and backend. From these Gaussian components, noise-free templates were 
created for each observing frequency with the phase of the fundamental component in the 
frequency domain rotated to zero. 

For each pulsar, the TOAs and ephemeris files were used in TEMPO to carry out a standard 
pulsar timing analysis (see Lorimer & Kramer 2005 and references therein). The Solar System 
model DE405 (Standish 1998) and the time standard TT(BIPM) were used in TEMPO. Time 
delays due to using two different backends were taken into account by placing jumps around 
the Spigot data, allowing a possible time offset due to each individual instrument. The result 
was a model for the behavior of the pulsar. The goodness of the fit was quantified by the 
root-mean-square of the difference between the model and the TOAs. The timing- derived 
parameters for each pulsar are listed in Tables 2 and 3 with doubled la errors listed to 
provide a conservative estimate of the 68% confidence limit. The root-mean-square values 
are on the order of a few milliperiods. Since no systematic trends are seen in the residuals, 
this is interpreted as an underestimation of the individual TOA errors. As is a common 
practice in pulsar timing, the TOA uncertainties are scaled by a constant factor, EFAC, to 
bring the reduced of each timing model to unity. 



^http: / /psrchive. sourceforge.net / current 
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4. 



Other Measured Parameters 



4.1. 



820-MHz Fluxes 



Fluxes were measured using the calibration routine associated with PSRCHIVE with the 
on-line folded GUPPl data. On- and off-source scans of the extra-galactic radio source 3C190 
were used for the flux calibration. Before each observation, a one-minute calibration scan 
was taken with a 25-Hz noise diode at the receiver. RFI was removed from the pulsar and 
calibration scan data using the psrzap utility from the PSRCHIVE package. 

PSRCHIVE calibration was applied to the data using on- and off-fiux scans, calibration 
scan, and pulsar data that allowed slightly different levels in the two polarizations (Hotan 
et al. 2004). These files were used to calculate an absolute flux conversion factor which was 
applied to the pulsar data. Lastly, a flux measurement was taken from the calibrated data. 



For pulsars with enough linearly polarized flux, rotation measures (RMs) were calcu- 
lated. This includes 10 of the 13 pulsars in this sample. PSRCHIVE was used with on-line 
folded GUPPI data for producing the RMs. The method used here was similar to that 
presented in the recent work of Yan et al. (2011). 

After the data were calibrated, we tested many different RM values from —1000 to 1000 
rad m~^ to find the one which provided the most polarized flux. An example plot showing the 
results of this procedure for PSR J194H-0121 can be seen in Figured A Gaussian was fit to 
the polarized flux versus RM and the peak of the Gaussian was taken as the best RM value. 
The RM values quoted for each pulsar are average values over multiple observations with the 
standard deviation as the uncertainty. This was done in case the ionosphere contaminated 
any RM measurements (see, e.g., Yan et al. 2011). 



Pulse widths at 50% intensity, 14^50, were derived either by fitting a Gaussian for single 
peak profiles or by measuring the width across the outer edge of the profile at 50% of the 
peak flux for multiple component profiles. This was done for each epoch. The W^o values 
for each epoch were averaged and the standard deviation was taken as the uncertainty. 



4.2. 



Rotation Measure 



4.3. 



Pulse Widths at 820 MHz 
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5. Discussion 

Results from pulsar timing and the analysis described in Section [Hare located in Tables 
2 and 3 for each individual pulsar and the location of each new pulsar on the P — P diagram 
is given in Figure HI The DM-derived distances are from the NE2001 electron density model 
(Cordes & Lazio 2002). The distances have uncertainties of approximately 25%. The surface 
magnetic field strength, characteristic age, and spin-down luminosity are calculated using 
standard formulae (see, e.g. Lorimer & Kramer 2005). Pulse profiles for each pulsar can 
be seen in Figures 6 and 7. Some of the individual pulsars are discussed in the following 
subsections. 



5.1. PSR J1327-0755 

At the time of discovery, the narrow 350-MHz pulse width of 0.25 ms, period of 2.68 ms, 
and location well off the Galactic plane (6 = 54.3°) of PSR J1327— 0755 were encouraging 
signs that this pulsar may be a good addition to pulsar timing arrays for the search for 
gravitational waves (e.g., Foster & Backer 1990). However, analysis at 820 MHz indicated 
that PSR J1327— 0755 may not meet the standards needed for a pulsar timing arrays pulsar 
due to very large variations in its flux from interstellar scintillation. While ~95% of GUPPI 
observations resulted in a detection, only ~35% of Spigot observations resulted in a detection. 
The pulsar was typically detected in only ~25% of the 200-MHz GUPPI bandpass, due to 
interstellar scintillation. The orbit has a period of 8.44 days and an eccentricity so low it 
has not yet been detected in our timing: e = 9(6) x 10^^. The projected semi-major axis of 
6.65 It-s implies a minimum companion mass of 0.22 M©, assuming a pulsar mass of 1.35 
Mq. Given the pulsar's location in the P — P diagram, the companion is very likely to be a 
white dwarf. 

A figure of merit for the strong equivalence principle (SEP; Damour & Schaefer 1991) 
test is the value of P^/e (Lorimer & Freire 2005). This value for PSR J1327-0755 is 7.91 x 
10^ days^. This values is within the range of values for other pulsars used for SEP tests and 
would be a good addition to arrays of pulsar used for SEP tests (Gonzalez et al. 2011). 

The composite proper motion for PSR J1327— 0755 is 99±23 mas yr^^. This proper 
motion corresponds to a transverse velocity of 800 ± 200 km s~^ for the predicted NE2001 
electron distance of 1.7 kpc. Only two recycled pulsars are known with larger values for 
proper motion and they are PSR J0437-4715 (Deller et al. 2009) and PSR J1231-1411 
(Ransom et al. 2011). Both of these pulsars are located within 0.5 kpc of the Earth. The 
large proper motion most likely implies that the NE2001 model does not accurately predict 
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distances off of the Galactic plane along this line of sight, or there are additional observational 
parameters that are not accounted for in PSR J1327— 0755's model. 

5.2. PSR J1623-0841 

PSR J1623— 0841 was discovered only in the single-pulse search with a period of 503 
ms and was labeled a rotating radio transient (RRAT; McLaughlin et al. 2006). Later 
observations at 350 MHz with Spigot showed single pulses, but no detections were made 
at 820 MHz. Gridding observations were unable to localize the position better than the 
discovery position because of its transient nature. The final TEMPO position was 20' away 
from the discovery position (outside the GBT's 820-MHz beam), and a timing solution was 
only attainable with the increased sensitivity of GUPPI and a dense set of observations at 
350 MHz. A full discussion of PSR J1623— 0841 will be presented elsewhere (Boyles et al., 
in preparation). 

5.3. PSR J1737-0814 

PSR J1737— 0814 is a 4.18 ms binary pulsar in a low-eccentric orbit (e = 0.000053) with 
an orbital period of 79.38 days and a minimum companion mass of 0.06 Mq, assuming a 
pulsar mass of 1.35 M0. Phinney (1992) predicted a relationship between orbital period and 
eccentricity of 

< >o-5= 1.5 X 10-^(P6/100 days) for Pb > Pent ~ 25 days, (2) 

based on fluctuation dissipation theory in the second star. The orbital eccentricity for 
PSR J1737— 0814 is about a factor of two lower than predicted by the model, but within the 
model's uncertainty. This indicates that the system may have followed the expected binary 
evolution for longer-period recycled pulsars. 

In principle, the low eccentricity should make this pulsar an excellent candidate for tests 
for violation of the Strong Equivalence Principle (Damour & Schaefer 1991) as well as the 
combined Lorentz-invariance and momentum-conservation strong-field Parameterized Post- 
Newtonian parameter (Bell & Damour 1996). However, the companion mass predicted by 
the core-mass-orbital-pcriod relationship (Tauris & Savonijc 1999) is approximately 0.33 Mq, 
meaning that, unless the system inclination angle is about 15° or less (which has an a priori 
likelihood of less than 4%), the system may have had a different evolutionary history. Since 
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it is best to use a set of sources that have followed the same evolutionary path for these tests 
(Wex 2000, Gonzalez et al. 2011), PSR J1737— 0811 may not be suitable for inclusion after 
all. We note that PSR J1327— 0755 has a minimum companion mass very close to the value 
predicted by the (Tauris & Savonije 1999) relation, therefore it is a pulsar that should be 
included in future editions of the ensemble tests for these parameters that would indicate a 
departure from General Relativity. 



5.4. PSR J1941+0121 

PSR J194H-0121 has the shortest period (217 ms) of the non-recycled pulsars found 
in this survey and about 85% of its total flux is linearly polarized. This high percentage 
of linearly polarized flux, along with its relatively wide pulse width (e = 0.093, where e is 
the pulse width divided by pulse period), allowed us to produce a model of the pulsar's 
emission using the PSRCHIVE software package. A complex-value Rotating Vector Model 
(RVM; Radhakrishnan & Cooke 1969) was fit to the Stokes Q and Stokes U profiles, treating 
them as real and imaginary numbers instead of the single parameter value of position angle. 
Figure |5] shows the best fit of the model and Table 1 gives the model fit parameters. The 
model gives an impact parameter (i.e. angle between magnetic field axis and line of sight, 
(3) of 8° ± 4° and a best-fit inclination angle (i.e. angle between rotation axis and magnetic 
field axis, a) of 138° ± 32° These parameter values indicate an inner line of sight for this 
pulsar. The model has a of 23.2 with 31 degrees of freedom (on-pulse phase bins minus 
number of fitted parameters), giving a reduced of 0.74. 



5.5. PSR J2222-0137 

PSR J2222— 0137 is in a low-eccentricity orbit (e = 0.00038) with an orbital period of 
2.4 days and a projected semi-major axis of 10.8 It-s. The orbital parameters imply a mass 
function of 0.229 M©. Assuming a pulsar mass of 1.35 Mq results in a minimum companion 
mass of 1.11 Mq. The spin period of 33 ms, the binary parameters, and the lack of a large P 
indicate that this pulsar is a partially recycled pulsar. Based on these parameters and what 
is known about other binary pulsars, this object is most likely a member of the so-called 
"intermediate-mass binary pulsars" (e.g. Camilo et al. 2001) in which the companion star is 
a CO white dwarf. It is also possible, though less likely, given the low orbital eccentricity, 
that PSR J2222— 0137 is a member of a double neutron star binary system. About a dozen of 
the early follow-up observations of PSR J2222— 0137 have been searched for a second pulsar 
without success with a minimum detectable flux density of 11 /iJy at 1400 MHz. 
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The presence of a second period derivative, suggests that there are unmodeled effects in 
the residuals. An unmodeled Shapiro delay signal (as the pulsar's emission passes through 
the gravitational potential of the companion) and no measured proper motion could be 
contributing to the second period derivative. We are not yet able to fit for this delay, but 
with a minimum companion mass of 1.11 M0 and narrow pulse duty cycle of 0.017, it 
should be measurable, leading to the companion mass and orbital inclination. Due to the 
possible effects of unmodeled variables and a reduced of 13.3, no EFAC is included for 
PSR J2222-0137's timing model, given in Table 3. An updated Shapiro delay and VLBA 
proper motion measurement, along with an updated timing solution, will be reported in a 
future publication. 

Using the DM value of 3.28 pc cm~^, the NE2001 electron density model predicts a 
distance of ~300 pc to PSR J2222-0137 (Cordes & Lazio 2002). If correct, PSR J2222-0137 
would be the second closest binary pulsar system after PSR J0437— 4715, which has a DM 
of 2.64 pc cm~^ and a distance of 156 pc (DcUcr ct al. 2009). For pulsars located this close 
to us, a large window of multi- wavelength observations opens up. These will include VLBA 
observations to measure proper motion and parallax, optical observations to search for an 
optical companion, X-ray observations to look for blackbody emission from the neutron star, 
and gamma-ray observations to look for higher magnetosphere emission. These observations 
will be reported in a future publication. 



6. Conclusions and Future Work 

The 350-MHz drift-scan survey revealed many interesting pulsars with a range of phys- 
ical questions that they can address. Along with PSR J2222— 0137, featured in this work, 
PSR J1023+0038 (Archibald et al. 2009), PSR J2256-1024 (Stairs et al. in prep), and 
PSR J0348-I-0438 (Paper II) have been studied in greater detail. Due to the success of this 
survey, a new survey named the Green Bank North Celestial Cap (GBNCC) survey is un- 
derway (Stovall et al. in prep). In this survey we are taking advantage of the newest pulsar 
instrument on the GBT, the GUPPI backend, which allows 100-MHz of bandwidth at 350 
MHz. This survey will be 30—40% more sensitive than the survey described in this paper 
and will be a powerful probe of the millisecond pulsar population in the northern sky. The 
GBNCC survey will eventually cover the whole entire northern sky visible by the GBT. 
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Fig. 1. — The sky coverage of the survey shown in Right Ascension and Dechnation (J2000 
coordinates). The Pulsar Search Collaboratory portion is shown in dark gray and the red 
dotted line shows the Galactic plane. New discoveries presented in this paper are shown as 
filled red triangles, while all other GBT drift-scan discoveries are shown as filled red stars. 



Parameter 


Angle 


Error 




(degrees) 


(degrees) 




-79 


3 




8 


4 


a 


138 


32 




351.7 


0.8 



Table 1: Geometric beam model fit for PSR J1941+0121 using a complex-number RVM, 
where 0o is the pulse phase of magnetic meridian, i/jq is position angle at (po, a is the 
colatitude of magnetic axis, and f3 is the impact parameter. 
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2 Pulses of Best Profile 



Candidate: ACCEUCand_1 
Telescope: GBT 

Epoch,--- = 54267.52159256296 




Epoch| 
T, 



'topo 
bory 



= 54267.52415597626 



Data Avg 
Data StdDev 
Profile Bins 
Profile Avg 
Profile StdDev 



0.00016384 

942080 

1.7Q3e+05 

683.3 

50 

3.208e+09 
9.379e+04 



Search Information 
RAj2(,„(, = 22:21:53.5575 DECj^ogo = -01:31:57.8026 

Best Fit Parameters 
Reduced = 7.464 P(Noise) < 4.99e-50 (•>='\A.Ba) 
Dispersion Measure (DM) — 3.290 
Ptopo ("^s) = 32.806019(95) 
P'lopo = 0.0(4.8)x10-= 

P"lopo (s/s^) = 0.0(2.0)x10 

Binary Parameters 
Pprb (s) = N/A e = N/A 

a,sin(i)/c (s) = N/A u (rod) = N/A 

T„„: = N/A 



P'Z. (s/s). 



bary 
bory 



32.809017(95) 
0.0(4.8)x10"^ 
(s/s^) = 0.0(2.0)x10-'° 




0.5 1 1.5 
Phase 

GBT350drift_54267_2221 



6 4 2 
Reduced 
-0131_DS2.fil 



DM 



0.01 -0.01 

Period - 32.80629777 (ms) 



4-Dec-2008 07:53 



Fig. 2. — Example of the search output diagnostic plot for new discovery PSR J2222— 0137. 
Left-hand side shows a pulse profile at the top and a plot of intensity versus time on the 
bottom. Along the plot of intensity versus time is a plot of accumulated versus time. 
The middle top plot shows observing frequency versus pulse phase and the bottom shows 
versus DM. The right hand-side (from top to bottom) shows versus period, versus P, 
and a intensity plot of period-P space. 
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1000 



tooo 



Rotation Measure [rad m^) 



Fig. 3. — Example fit for RM for PSR J194H-0121. The data points are square boxes with 
error bars and the best-fit RM is shown by a vertical dotted line. 
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Fig. 4. — A P — P diagram showing the new discoveries reported in this paper and the 
companion paper Paper II as black diamonds. The lines on the plot show constant spin-down 
luminosity (dotted line), constant magnetic field (dash line), and constant characteristic age 
(dot-dash line). The pulsar population is taken from the ATNF pulsar catalog (Manchester 
et al. 2005, ,http://www. atnf.csiro.au/research/pulsar/psrcat/ ). 
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Fig. 5. — RVM model fit for tlie radio emission of PSR J1941+0121. The top plot shows 
parallactic angle versus pulse phase, while the bottom plot shows integrated pulse flux density 
versus pulse phase for total intensity (black line), linearly polarized flux (red line), and 
circularly polarized flux (blue line). 



Parameter 



PSR J1555-0515 PSR J1612+2008 



Timing Parameters 



Right Ascension (J2000) 


15:55:40.097(16) 


16:12:23.432(3) 


Declination (J2000) 


-05:15:57.4(5) 


20:08:18.33(4) 


Spin Period (s) 


0.97540989108(5) 


0.4266459810253(12) 


Period Derivative (s s~^) 


2.054(2) X 10-^^ 


3.727(15) X 10-^^ 


Dispersion Measure (pc cm"^) 


23.46(3) 


19.544(10) 


Reference Epoch (MJD) 


54985 


54985 


Span of Timing Data 


54635-55335 


54635-55335 


Number of TOAs 


36 


41 


Weighted RMS Residual {/is) 


524 


164 


EFAC 


3.22 


1.72 


Derived Parameters 


Galactic Longitude (degrees) 


3.97 


35.52 


Galactic Latitude (degrees) 


34.97 


43.74 


Distance (kpc) 


1.2 


1.3 


Distance From Plane (kpc) 


0.69 


0.90 


Surface Magnetic Field (lO"*^^ Gauss) 


1.4 


0.13 


Spin-down Luminosity (10^^ ergs s~^) 


8.8 


1.9 


Characteristic Age (Myr) 


7.4 


180 


820 MHz Flux Density (mJy) 


1.2(2) 


0.83(11) 


W50 820 MHz (ms) 


12.9(3) 


3.29(10) 


Rotation Measure (rad m~^) 


1.3(3.0) 


22(3) 



Table 2: Timing and derived parameters for newly discovered isolated pulsars. 
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Parameter 


PSR J1623-0841 


PSR J1633-2009 




Timing Parameters 


Right Ascension (J2000) 


16:23:42.701(19) 


16:33:55.30(6) 


Declination (J2000) 


-08:41:36.4(9) 


-20:10:09(5) 


Spin Period (s) 


0.503014997755(14) 


0.93555704483(6) 


Period Derivative (s s~^) 


1.958(3) X 10-^5 


1.070(3) X 10-^5 


Dispersion Measure (pc cm"^) 


60.42(4) 


48.19(6) 


Reference Epoch (MJD) 


55079 


54993 


Span of Timing Data 


54635-55522 


54651-55335 


Number of TOAs 


36 


29 


Weighted RMS Residual (/xs) 


670 


633 


EFAC 


2.17 


1.62 


Derived Parameters 


Galactic Longitude (degrees) 


5.77 


357.63 


Galactic Latitude (degrees) 


27.37 


18.33 


Distance (kpc) 


3.3 


1.6 


Distance From Plane (kpc) 


1.5 


0.5 


Surface Magnetic Field (lO"*^^ Gauss) 


1.0 


1.0 


Spin-down Luminosity (10^^ ergs s~^) 


63 


5.0 


Characteristic Age (Myr) 


4.1 


14 


820 MHz Flux Density (mJy) 


0.42(13) 


1.35(4) 


W50 820 MHz (ms) 


11.2(3) 


32.2(1.0) 


Rotation Measure" (rad m~^) 




-0.1(1.2) 



" PSR J1623— 0841 does not have enough hnearly polarized flux for a RM measurement. 



Table 2: Timing and derived parameters for newly discovered isolated pulsars. 



-23- 



Parameter PSR J1735-0243 PSR J1903-0848 



Timing Parameters 



Right Ascension (J2000) 


17:35:48.1(2) 


19:03:11.271(18) 


Declination (J2000) 


-02:43:48(5) 


-08:48:57.4(8) 


Spin Period (s) 


0.7828869765(5) 


0.88732464056(4) 


Period Derivative (s s~^) 


6.32(12) X 10-^^ 


1.330(2) X 10-^5 


Dispersion Measure (pc cm~^) 


55.4(5) 


66.99(4) 


Reference Epoch (MJD) 


54827 


54987 


Span of Timing Data 


54352-55302 


54634-55339 


Number of TOAs 


39 


51 


Weighted RMS Residual (/xs) 


6532 


824 


EFAC 


3.70 


3.32 


Derived Parameters 


Galactic Longitude (degrees) 


21.56 


26.38 


Galactic Latitude (degrees) 


15.45 


-6.60 


Distance (kpc) 


1.9 


2.0 


Distance From Plane (kpc) 


0.51 


0.23 


Surface Magnetic Field (lO"*^^ Gauss) 


0.68 


1.1 


Spin-down Luminosity (10^^ ergs s~^) 


4.7 


7.2 


Characteristic Age (Myr) 


21 


11 


820 MHz Flux Density (mJy) 


1.3(3) 


1.4(2) 


W50 820 MHz (ms) 


102(2) 


12.3(2) 


Rotation Measure (rad m~^) 


28.3(1.9) 


4.4(9) 



Table 2: Timing and derived parameters for newly discovered isolated pulsars. 



Parameter 



PSR J1941+0121 PSR J2012-2029 



Timing Parameters 



Right Ascension (J2000) 


19:41:16.039(16) 


20:12:46.6(3) 


Declination (J2000) 


01:21:39.5(5) 


-20:29:31(22) 


Spin Period (s) 


0.217317451828(10) 


0.54400187061(3) 


Period Derivative (s s~^) 


1.913(5) X 10-^^ 


5.475(18) X 10-^6 


Dispersion Measure (pc cm"^) 


52.7(7) 


37.67(8) 


Reference Epoch (MJD) 


55026 


54987 


Span of Timing Data 


54712-55339 


54635-55339 


Number of TOAs 


33 


33 


Weighted RMS Residual (/xs) 


555 


538 


EFAC 


1.98 


1.59 


Derived Parameters 


Galactic Longitude (degrees) 


39.91 


22.41 


Galactic Latitude (degrees) 


-10.43 


-26.70 


Distance (kpc) 


2.2 


1.4 


Distance From Plane (kpc) 


0.40 


0.63 


Surface Magnetic Field (lO"*^^ Gauss) 


0.20 


0.55 


Spin-down Luminosity (10^^ ergs s~^) 


72 


14 


Characteristic Age (Myr) 


18 


16 


820 MHz Flux Density (mJy) 


1.7(5) 


1.00(12) 


W50 820 MHz (ms) 


20.2(1.1) 


24.2(7) 


Rotation Measure" (rad m~^) 


-92.0(3) 





" PSR J2012— 2029 does not have enough hnearly polarized flux for a RM measurement. 



Table 2: Timing and derived parameters for newly discovered isolated pulsars. 
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Parameter PSR J2033-1938 PSR J2111+2106 



Timing Parameters 



Right Ascension (J2000) 


20:33:55.4(2) 


21:11:33.13(6) 


Declination (J2000) 


-19:38:59(12) 


21:06:07.0(1.4) 


Spin Period (s) 


1.28171902379(12) 


3.9538529596(6) 


Period Derivative (s s~^) 


4.55(6) X 10-^^ 


3.24(2) X 10-^5 


Dispersion Measure (pc cm~^) 


23.47(9) 


59.74(14) 


Reference Epoch (MJD) 


54987 


54987 


Span of Timing Data 


54635-55339 


54635-55339 


Number of TOAs 


37 


35 


Weighted RMS Residual (/^s) 


606 


2106 


EFAC 


2.05 


2.33 


Derived Parameters 


Galactic Longitude (degrees) 


25.32 


69.40 


Galactic Latitude (degrees) 


-31.04 


-18.20 


Distance (kpc) 


1.00 


3.8 


Distance From Plane (kpc) 


0.51 


1.2 


Surface Magnetic Field (lO"*^^ Gauss) 


0.76 


3.6 


Spin- down Luminosity (10'^^ ergs s~^) 


0.85 


0.20 


Characteristic Age (Myr) 


44 


19 


820 MHz Flux Density (mJy) 


1.13(15) 


1.08(10) 


W50 820 MHz (ms) 


32.4(3) 


45(2) 


Rotation Measure (rad m~^) 


-17.7(5) 


-75.3(8) 



Table 2: Timing and derived parameters for newly discovered isolated pulsars. 
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Parameter 


PSR J1327-0755" 


PSR J1737-0811 




Timing Parameters 


Right Ascension (J2000) 


13:27:57.5880(12) 


17-37-47 11235fl4l 

-LI •tJl •^l - J — I- \ J_ ^ 1 


Declination (J2000) 


-07:55:29.80(4) 


-08:11:08.887(5) 


RA Proper Motion 


27(18) 




DEC Proper Motion 


95(48) 




Spin Period (s) 0.0026779231971205(2) 


0.0041750173128551(8) 


Period Derivative'' (s s ^) 


1.773(18) X 10-20 


7.93(8) X 10-21 


Dispersion Measure (pc cm~'^) 


27.91215(6) 


55.311(3) 


Reference Epoch (MJD) 


55131 


54987 


Span of Timing Data 


54649-55613 


54635-55339 


Number of TOAs= 


188 


225 


Weighted RMS Residual (iis) 


5.5 


14 


EFAC 


1.33 


1.08 


Binary Parameters 


Binary Model 


ELLl 


BT 


Orbital Period (days) 


8.439086019(12) 


79.517379(2) 


Projected Semi-major Axis (It-s) 


6.645774(2) 


9.332791(3) 


Time of Ascending Node (MJD) 


54717.3830798(7) 




First Laplacc-Lagrangc Parameter 


6(5) X lO-'^ 




Second Laplace-Lagrange Parameter 


6(7) X lO-'^ 




Epoch Of Periastron (MJD) 54718.50769778(9 x 10"^) 


54696.879781933(2x10-1) 


Orbital Eccentricity 


9(6) X 10-^ 


5.38(8) X 10-5 


Longitude of Periastron (degrees) 


48(40) 


49.8(9) 


Mass Function (Mg) 


0.004425158(4) 


0.0001380362(2) 


Dcri\"(Hl Puraiuct c>rs 


Minimum Companion Mass'' ( Mq) 


0.2219585(9) 


0.0662357(3) 


Galactic Longitude (degrees) 


318.39 


16.86 


Galactic Latitude (degrees) 


53.85 


12.34 


Distance (kpc) 


1.7 


1.7 


Distance From Plane (kpc) 


1.0 


0.36 


Surface Magnetic Field'' (10^ Gauss) 


0.70 


1.8 


Spin-down Luminosity'' (10^^ ergs s~^) 


360 


430 


Characteristic Age'' (Gyr) 


2.4 


8.3 


820 MHz Flux Density (mJy) 


0.9(2) 


2.7(2) 


W50 820 MHz (ms) 


0.131(5) 


0.52(6) 


Rotation Measure^ (rad m~^) 




71.4(1.4) 



Epoch Of Periastron, Orbital Eccentricity, and Longitude of Periastron are derived 



values for PSR 1327-0755. 

'' Values are not corrected for Shklovskii effect (Shklovskii 1970). 

Three TOAs produced per observation. 

Assumes a pulsar mass of 1.35 Mq. 
* PSR J1327— 0755 does not have enough linearly polarized flux for a RM measurement. 

Table 3: Timing, binary, and derived parameters for newly discovered binary pulsars. 
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Parameter 


PSR J2222-0137 




Timing Parameters 


Right Ascension (J2000) 


22:22:05.96345(8) 


Declination (J2000) 


-01:37:15.627(3) 


Spin Period (s) 0.032817859050650(5) 


Period Derivative (s s~^) 


4.74(3) X 10-20 


Second Period Derivative (s s~^) 


1.61(8)xl0-27 


Dispersion Measure (pc cm~^) 


3.27511(10) 


Reference Epoch (MJD) 


55285 


Span of Timing Data 


54985-55585 


Number of TOAs" 


219 


Weighted RMS Residual {/is) 


8.2 


Binary Parameters 


Binary Model 


BT 


Orbital Period (days) 


2.4457599960(5) 


Projected Semi-major Axis (It-s) 


10.8480469(8) 


Orbital Eccentricity 


0.00038454(14) 


Longitude of Periastron (degrees) 


119.863(14) 


Time of periastron (MJD) 55284.775883844(1.0 x 10"^) 


Mass Function ( M,,,) 


0.22911166(5) 


Derived Parameters 


Minimum Companion Mass'' ( Mq) 


1.1169752(3) 


Galactic Longitude (degrees) 


62.02 


Galactic Latitude (degrees) 


-46.08 


Distance (kpc) 


0.31 


Distance From Plane (kpc) 


0.22 


Surface Magnetic Field (10^ Gauss) 


1.3 


Spin-down Luminosity (10^^ ergs s~^) 


5.3 


Characteristic Age (Gyr) 


11 


820 MHz Flux Density (mJy) 


2.6(5) 


Wso 820 MHz (ms) 


0.570(5) 


Rotation Measure (rad m^^) 


1.8(6) 


° Three TOAs produced per observation. 




^ Assumes a pulsar mass of 1.35 M©. 





Table 3: Timing, binary, and derived parameters for newly discovered PSR J2222— 0137. 
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J1555-0515 
975 ms 



J1612+2008 
427 ms 



J1623-0841 
503 ms 



J1633-2009 
936 ms 



J1735-0243 
783 ms 



J1903-0848 
887 ms 



J8012-2029 
544 ms 



J2033-1938 
1282 ms 



J2111+2108 
3954 ms 



Fig. 6. — Profiles at 820-MHz for nine newly discovered pulsars created from multi-epoch 
data using the GUPPI backend. Spin periods are listed. 
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Fig. 7. — Profiles at 820-MHz for the remaining four newly discovered pulsars created from 
multi-epoch data using the GUPPI backend. The spin periods are hsted. 



